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1.  INTRODUCTION 

In  the  theory  of  heterogeneous  catalysis,  it  is  generally  assumed 
that  the  d electrons  of  transition  metal  catalysts  are  dominant  in  the 
desired  chemical  reaction  taking  place  near  the  surface.^  However,  a 
microscopic  mechanism  describing  how  the  d electrons  in  the  surface 
states  enter  into  the  chemical  reaction  >-has  not  previously  been 
proposed.  Also,  the  f electrons  are  assumed  responsible  for  the 
successful  catalytic  activity  of  compounds  composed  of  rare-earth  ions; 
a successful  theory  is  lacking. 

We  have  proposed  a concept  describing  how  a rare-earth  ion  in  the 

surface  of  a metal  might  enter  into  a reaction  of  an  atom  near  the 

surface  and  increase  the  rate  of  ionization  of  the  atom  under  certain 

favorable  conditions.  The  concept  is  based  upon  an  extension  of  a 

previous  theory  of  the  crystal  field  parameters,  which  explained 

the  splitting  of  energy  levels  of  rare-earth  ions  in  bulk  rare-earth 

monochalcogenides  and  monopnictides . That  theory  was  based  upon  the  use 

of  the  Thomas-Fermi  model  of  electrons  in  a metal  and  a sum  over  the 

various  constituent  ions  to  obtain  the  lattice  sum  parameters,  A , and 

nm 

consequently  the  B . Then  the  lattice  sum  was  formed  for  an  ion  in  the 
nm 

surface  of  the  metal,  and  the  B for  an  ion  in  the  surface  were 

nm 

obtained.  With  these  parameters,  the  energy  levels  and  corresponding 
wave  functions  were  obtained.  These  wave  functions  were  used  in  a 

semiclassical  approximation  to  calculate  the  charge  density  of  the 
rare-earth  ion  and,  consequently,  the  electric  potential  near  the 
surface  of  the  metal.  This  potential  was  shown  to  depend  on  the  energy 

^J.  C.  Slater  and  K.  H.  Johnson,  Quantum  Chemistry  and  Catalysis , 
Physics  Today  (October  1974),  34. 

^C.  A.  Morrison , N.  Karayianis , and  D.  E.  Wortman,  A Possible  Use  of 
the  Surface  States  of  Transition  and  Rare-Earth  Metal  Ions  in  the  Theory 
of  Catalysis , Harry  Diamond  Laboratories  TR-1752  (April  1976). 

^C.  A.  Morrison,  Phys.  Lett.,  51A  (1975) , 49. 

5 


I 


[ 


t 


level  of  the  rare-earth  ion  in  the  surface  and  also  to  depend  on  the 

crystallographic  plane  of  the  surface.  This  potential  along  with  a 

calculation  using  standard  metal-ion  interactions  was  proposed  to 
illustrate  the  utility  of  the  theory.  For  those  calculations,  we  chose 
two  papers  by  Gadzuk,**'^  in  which  the  ionization  of  alkali  earths  by  W 
or  Pt  was  calculated  by  using  perturbation  theory  and  metal-ion 
interactions  obtained  by  using  classical  electrical  image  theory. 

In  this  report,  we  investigate  Gadzuk's  theory  for  his  examples  and 
correct  some  apparently  typographical  errors  in  his  papers.  Then  we 
extend  his  results  to  include  the  interactions  created  by  the  ionic 

nature  of  the  atom  in  the  surface  of  the  metal.  Finally,  this 

calculation  is  done  for  the  metal  PrSb  and  the  added  interaction  created 
because  the  rare-earth  ion  in  the  surface  is  in  one  of  the  low-lying 
surface  energy  levels. 

2.  GADZUK’S  THEORY  OF  SURFACE  INTERACTIONS 

Gadzuk’'"^'  gave  a theory  of  ad-atom-metal  interactions,  in  which  the 
ad-atom  is  assumed  to  interact  with  the  metal  according  to  classical 
image  theory.  Details  of  the  calculation  of  the  energy  shift  are  given 
in  appendix  A.  Here  we  state  the  underlying  physical  considerations  and 
give  the  results  for  Cs  adsorbed  on  W and  K adsorbed  on  Pt,  which  Gadzuk 
considered. 

2 . 1 Interactions  of  Atom  and  Solid 

Consider  the  physical  configuration  and  coordinate  system 
shown  in  figure  1.  The  nucleus  of  the  ad-atom  is  at  the  point  (0,0,0); 

“*J.  W.  Gadzuk,  Surface  Science,  £ (1967),  133. 

^J.  W.  Gadzuk,  Surface  Science,  ^ (1967),  159. 
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Figure  1.  Physical  configuration  and  coordinate  system  for  alkali 
earth  adsorption. 


the  electron  is  at  r = (x,y,z).  The  half  space  z < -s  is  assumed  to  be 
occupied  by  metal.  The  interaction  of  the  atom  with  the  metal  is 
represented  by  classical  image  theory  which,  for  the  electron  on  the 
alkali  earth,  is  given  by  (app  A) 


V. 

1 


(1) 


v^ere 


V. 

1 

is 

the 

interaction 

energy , 

q^/4di 

is 

the 

interaction 

of  the  electron  with 

its 

own 

image , 

q^/R 

is 

the 

interaction 

between  the  electron 

and 

the 

ion  image, 

q2/r 

is 

the 

interaction 

of  the  electron  with 

the 

ion 

core . 

This  last  term  is  present  even  in  the  cibsence  of  the  metal.  Gadzuk  thus 
chooses  the  single  electron  Hamiltonian,  H,  to  be 


J 


tl  0 

— V'^-V,z<s  -s, 

2m  o c 


z>s^-s 


Where  s is  chosen  so  that 
c 


2 2 

- 2_  + ^ 2 = s _ g 

4di  R c 


vdiere 


fi  is  Planck's  constant, 

is  the  electron  mass, 
is  the  Laplacian  operator, 

V is  the  sum  of  the  Fermi  energy,  E , and  the.  work  function,  <t, 
o F 

z is  the  coordinate  normal  to  the  metal  surface, 

s is  the  cutoff  distance, 

c 

s is  the  distance  of  the  nucleus  from  the  surface. 

If  we  assume  that  the  coordinates  x and  y of  the  electron  are  small 

compared  to  s (the  assumption  made  by  Gadzuk)  the  value  of  s^  given  by 
equation  (3)  is 


j jj(s  + A)(s  + 91)1^  - s - 6^|, 


where  \ = q^/4V  . If  V is  given  in  electron  volts,  then  X = 3.6/V  in 
o o o 

angstrom  units.  The  variation  of  s^  with  s is  shown  in  figure  2 for  Pt 
and  figure  3 for  W.  (The  curves  labeled  by  ionic  charge  Qj  = 0 are 
pertinent  here.) 
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Figure  2.  Cutoff  parameter  for  platinum. 


2.2  Energy  Shift 


In  his  theory  of  ionization  of  alkali  earths , Gadzuk  chose  as 
wave  functions  of  the  electron  the  2s  hydrogenic  type  wave  functions 

i|i  = N(1  - ar)  e (5) 


with  01  = 0.99  A 


for  Cs  and  a = 1.16  A for  K.  These  values  of 


a were  chosen  so  that  the  hydrogenic  wave  functions  closely 

approximated  the  exact  self-consistent  calculations  of  Herman  and 
Skillman. ^ The  energy  shift  of  the  electron  on  the  alkali  atom  is 
obtained  by  using  first-order  perturbation  theory.  That  is, 

®F.  Herman  and  S.  Skillman,  Atomic  Structure  Calculations,  Prentice 
Hall,  Inc.,  Englewood  Cliffs,  NJ  (1963). 


where 


VERSUS  S 
Cs  ON  W 
E,  +®  = 10,3 
a = 3.16469  k 
!100I  PUNE 


DISTANCE  OF  NUCLEUS  FROM  SURFACE,  s (A) 


Figure  3.  Cutoff  parameter  for  tungsten. 
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AE^  is  the  energy  shift. 


* stands  for  complex  conjugation, 


covers  the  region  -®<x,y<®,s  - s<z<°o, 

c 


and  where 


V.  = — a! — aL_ 

1 4(s  + z)  2s  + z 


The  assumption  R 2s  + z has  been  made  in  equation  (7),  By  using  the 
wave  function  given  by  equation  (5),  the  integrals  of  equation  (6)  may 
be  evaluated.  The  result  is 


2 ( 

AE^  = 3—“  |4g  (2as)  ^Ei(4as)  - Ei^2a(s  + s^)j| 

-4g  (-2as)Ei  (-4as)  - g (as)  ^Ei  (2us ) - Ei(2as^)J 

+g  (-as ) Ei(-2as ) + |"  - s^,)  + 25(as)^ 

-ISa^ss^  + 3 (asg)  ^Jexp  [-2a  (s  - Sj,)]|/F(s,  s^) 


where 


-(i 


+ X - + 2x^1 e 


F(s,  s^) 


2 - 1 1 + ja  (s  - s^)  + a^(s  - s^,)  ^ + a^(s 


• exp  [-2a  (s  - s^)]  , 


Ei(z) 


=r4- 


n 

as  derived  from  Abramowitz  and  Stegun.  The  details  of  the  derivation 
are  given  in  appendix  A.  This  result,  equation  (8) , is  to  be  compared 
with  equation  (11)  in  Gadzuk's  paper,  vdiich  is  incorrect.  Gadzuk's 
result  should  also  be  multiplied  by  two  if  the  energy  shift  is  to  be 
given  in  electron  volts,  and  his  figure  4 should  have  the  ordinate 
doubled. 


**J.  W.  Gadzuk,  Surface  Science,  ^ (1967),  133. 

Abramowitz  and  I.  A.  Stegun,  Ed.,  Handbook  of  Mathematical 
Functions , National  Bureau  of  Standards , Gaithersburg , HD  (1964)  . 
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Figure  4.  Energy  shift  in  Gadzuk's  theory. 

The  results  given  in  equation  (8)  are  shown  in  figure  4,  which 
should  be  compared  to  Gadzuk's  figure  4.  These  results  are  quite 
similar  when  the  factor  of  two  is  taken  into  consideration.  Slight 
changes  in  the  Fermi  energy  or  work  function  can  account  for  the 
difference  in  the  results.  (Gadzuk  did  not  state  the  values  of  the 
Feirmi  energy  or  work  function  used  in  his  calculations.) 


The  energy  shift  given  above  is  that  predicted  on  the  basis  of  jjj 

classical  image  theory.  We  shall  now  consider  the  effects  of  several 
refinements  of  the  theory. 


3.  ENERGY  SHIFT  DUE  TO  IONIC  CHARGES  AT  SURFACE 
3 . 1 Potential 

If  we  assume  that  the  ion  cores  of  the  surface  atoms  carry  a 
residual  charge  which  is  neutralized  by  a continuous  background  of 
electrons,  then  the  results  in  section  2 are  significantly  altered.  In 
the  work  reported  here,  we  assume  that  only  these  ions  on  the  surface 
contribute  to  the  potential  outside  the  metal.  The  charge  distribution 
in  the  surface  (the  plane  z = -s)  is  assumed  to  be 


o(x,y)  =^Ql'5(x  - na)6(y  - 


where 


a is  the  surface  charge  density, 
x,y  are  coordinates  in  the  surface  plane. 


a is  the  lattice  constant. 


(5  is  the  Dirac  delta  function,  defined  by 


J 6 (x  - x^^f(x)dx  = 


for  body-centered  cubic  W and 


(x,y)  Ql|<5(x  - na)6(y  - ma) 

nm  ( 

+ 6[x  - (n  + |)a]6[y  - (m  t |]a]j - 


for  face-centered  cubic  Pt.  The  first  terms  involving  6 -functions 
represent  the  ionic  point  charges,  and  the  last  term  represents  the 


13 


continuous  electron  distribution.  The  neutrality  of  the  system  can  be 

demonstrated  by  integrating  either  distribution  (eg  (12)  or  (13))  over 

the  area  occupied  by  a single  cell.  The  potential,  V (z) , at  the 

P 

electron  on  the  alkali  atom  is  derived  in  appendix  B.  The  cutoff 
condition  given  by  equation  (3)  is  now  modified  to 


2 2 

-V  = + ^ + V (z)  , z = s - s, 

o 4dj  R p c 


(14) 


which  determines  s . Equation  (14)  has  been  used  to  determine  s as  a 
c c 

function  of  s for  several  values  of  (fig.  2,  3).  The  addition  of  the 

potential  created  by  the  surface  ions  tends  to  lower  the  barrier  near 

the  surface  of  the  metal  and  thereby  decreases  the  total  potential  to 

-V  at  greater  distances  from  the  surface, 
o 

3.2  Energy  Shift 

The  energy  shift  given  by  equation  (8)  is  considerably 

modified  by  the  inclusion  of  the  ions  in  the  surface  of  the  metal.  The 
surface  ions  have  a large  effect  on  the  cutoff  near  the  surface 

(fig.  2,  3) . 


A difficulty  exists  in  determining  the  value  of  the  ionic 
charge,  and,  for  the  purposes  here,  we  have  chosen  both  Pt  and  W to  have 
a charge  of  +e.  Much  larger  values  would  be  necessary  for  W if  the 
charge  were  chosen  proportional  to  the  valence  of  the  metal  ion. 

The  calculation  of  the  energy  shift  due  to  the  point-charge 
lattice  on  the  surface  is  straightforward  (app  C) . The  result, 
equation  (C-13) , must  be  added  to  the  image  contribution,  equation  (8) , 
to  obtain  the  total  energy  shift  (fig.  5)  for  K on  Pt  or  Cs  on  W. 
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OISTANCE  OF  NUCLEUS  FROM  SURFACE  (A| 

Figure  5.  Energy  shift  including  point  charge  potential. 


The  energy  shift  may  be  decomposed  into  three  parts: 

a.  The  image  shift  as  calculated  by  Gadzuk 

b.  The  change  in  Ae.  due  to  the  change  in  s , as  given  by 

i c 

equation  (14) 

c.  The  contribution  of  the  point-charge  potential 

We  have  plotted  these  components  separately  for  K adsorbed  on 

Pt,  The  image  shift  is  given  in  figure  4.  Figure  2 shows  the  change  in 

s for  several  values  of  Qi.  The  effect  of  this  change  on  the  energy 
c 
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shift  is  shown  in  figure  6.  Finally,  the  contribution  of  the  point 
charges  directly  is  shown  in  figure  7.  For  Cs  adsorbed  on  W,  we  have 
given  the  effect  of  the  point-charge  potential  in  changing  the  cutoff  in 
figure  3. 


When  considering  the  detailed  potential  near  the  surface  of  a 
metal,  we  must  allow  for  the  structure  of  the  surface  states;  in 
particular,  there  arise  contributions  to  the  potential  due  to  the 
asphericity  of  the  surface  states.  These  contributions  are  considered 
in  section  4 for  a rare-earth  metal.  It  is  this  contribution  that  may 
distinguish  catalytic  metals  from  ordinary  metals. 


9 10 
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SHIFT  DUE  TO  IMAGE  POTENTIAL 
Ef  + ct>  = IP  2 eV 


Figure  b.  Change  in  energy  shift  due  to  change  in  cutoff. 


We  now  consider  those  contributions  to  the  potential  near  the 
surface  of  a metal  which  depend  on  the  detailed  surface  states  of  the 
ions  at  the  surface.  For  a metal  composed  of  rare-earth  ions,  this 
contribution  may  be  done  simply,  since  the  f-electrons  may  be  regarded 
as  localized  about  particular  ionic  sites.  This  localization 
corresponds  to  the  tight-binding  approximation  in  band  theory.^  We 
shall  take  the  approach  of  crystal  field  theory.  That  is,  the  electric 
potential  near  a rare-earth  site  is  given;  the  energy  levels  and  wave 


R.  Harrison,  Pseudopotentials  in  the  Theory  of  Metals,  Benjamin, 

I New  York  (1966) . 

* 9 

B.  G.  Wybourne,  Spectroscopic  Properties  of  Rare  Earths, 
] Jntersciences  Publishers , New  York  (1965) . 

I 
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functions  of  the  rare-earth  ion  in  the  crystal  field  are  then  computed, 
and  from  them  the  multipolar  moments  of  the  rare  earth  can  be 
determined.  These  multipolar  moments,  qj^^»  depend  on  the  state  of  the 
rare-earth  ion.  The  potential  near  the  surface  may  then  be  represented 
in  terms  of  a lattice  sum  involving  the  ionic  multipolar  moments  over 
the  surface  of  the  metal.  Detailed  calculations  of  this  contribution  to 
the  potential  are  given  in  appendix  B. 

We  have  calculated  the  potential  and  energy  shift  for  Cs  and  K 
adsorbed  on  PrSb,  which  earlier  was  illustrated  for  this  type  of 
potential  near  the  surface.^  Multipolar  moments  of  the  Pr^"*”  ion  in 
various  states  of  the  ^H4  ground  multiplet  have  been  computed  (table  I) . 

TABLE  I.  MATRIX  ELEMENTS  . I r"c  IP.  >FOR  THE 

\ 1 ' nm'  1/  4 

MULTIPLET  OF  Pr^'*'  IN  THE  SURFACE  OF  PrSb 
(A  UNITS) ^ 


Energy 

r 

n,m 

2.0 

4,0 

4,4 

6,0 

6,4 

0 

1 

0.06770 

-0.02043 

-0.00544 

-0.03075 

0.00730 

49.5 

3 

0.05850 

-0.01065 

-0.00B94 

0.00037 

-0.00229 

150.0 

2 

0.03277 

0.00961 

-0.03633 

0.03056 

-0.03627 

252.0 

2 

0.03050 

0.01031 

0.03451 

0.03046 

0.03621 

415.0 

3 

-0.01659 

0.02445 

0.00661 

-0.01895 

0.00413 

592.0 

1 

-0.08486 

-0.01663 

-0.00096 

0.00390 

0.00064 

595.0 

1 

-0.08417 

-0.01673 

0.00220 

0.00380 

-0.00B76 

' for  Pr'* , 

th»  follot^ing  vslu9M  of 

hsv9  boon 

usod: 

<r’>  - 0.304194  (A)' 

<r*>-  0.13131}  (A)' 

<r‘>  - 0.343410  (A)^ 

A.  J.  Fre^mAft  *ruS  A.  B.  Vstson»  Phys . kev,,  127  (19f>2) , 20^9. 


^C.  A.  Morrison,  N.  Karaqianis , and  D.  E.  Wortman,  A Possible  Use  of 
the  Surface  States  of  Transition  and  Rare-Earth  Metal  Ions  in  the  Theory 
of  Catalysis,  Harry  Diamond  Laboratories  TR-1752  (April  1976). 
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These  have  been  used  with  the  results  of  appendix  C (eq  (C-14) ) to 
obtain  energy  shifts  due  to  the  surface  states.  These  additional  shifts 
are  shown  in  figure  8 for  K and  figure  9 for  Cs. 

Figure  10  shows  the  effect  of  the  surface  states  in  changing  the  s 

c 

for  the  ground  state  of  Pr.  The  condition  used  to  obtain  s^  is  similar 
to  equation  (14),  except  that  the  surface-state  potential,  V^,  is  added 
to  the  right-hand  side.  This  addition  results  in  a slight  change  in  s^ 
(fig.  10) . Figure  11  shows  the  total  energy  shift  for  K adsorbed  on  the 
[100]  surface  of  PrSb.  This  energy  shift  includes  all  the  effects 
considered  in  sections  2 and  3,  as  well  as  the  surface-state 
contributions  considered  here,  with  the  Pr  ions  in  their  ground  state. 
Figure  12  gives  a similar  plot  for  Cs  adsorbed  on  PrSb. 


Figure  8.  Energy  shift  in  [100]  plane  due  to  surface  states 
for  K on  PrSb. 


I'l 


NOTE  NUMBEfIS  IN  PARENTHESES  ARE 
ENERGY  VALUES  IN  cm  ' 


Figure  9.  Energy  shift  n [100]  plane  due  to  surface  states 
for  Cs  on  PrSb. 


PrSb 


Figure  10.  Cutoff  parameter  for  PrSb. 


This  section  completes  the  description  of  the  surface  potential  of 
a metal  according  to  Gadzuk's  theory  (with  the  additions  of  the 
point-charge  and  surface-state  potentials).  In  section  5,  we  point  out 
some  deficiencies  in  Gadzuk's  theory,  their  effects  on  our  calculations, 
and  remedies  for  them. 


5.  CRITICISM 


The  formalism  introduced  by  Gadzuk  for  calculating  the  energy 

shifts  of  alkali  earth  atoms  near  metals  is  unsatisfactory  in  many 

respects.  First,  the  is  not  a good  representation  of  the  imperfect 

screening  of  applied  electric  fields  by  the  metal.  Second,  Gadzuk 

neglected  some  important  contributions  to  the  energy  shift  which  would 

be  present  in  a consistent  treatment.  Third,  although  Gadzuk  claims 

that  his  procedure  is  first-order  perturbation  theory,  the  assumption 

that  the  wave  function  vanishes  for  z < s - s is  not  standard. 

c 

5.1  Image  Potential 


The  image  potential  given  by  Gadzuk  is  valid  actually  only  at 
large  distances  from  the  metal.  Other  attempts  have  been  made  to  cast 
the  image  potential  in  a form  which  is  approximately  correct  everywhere. 
One  of  these  is  ' 


V(r) 


4 (s  + s + z) 
c 


z > 0 , 


V , z < 0 , 

o 


(15) 


M.  Propst,  Phys.  Rev.,  129  (1962),  7. 
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in  vrtiich  the  is  explicit  in  the  potential  and  is  given  by 

s^  = qV4V^  . (16) 

Energy  shifts  may  be  calculated  by  using  equation  (15)  in  a manner 
similar  to  that  used  in  section  2.  However,  the  conceptual  difficulties 
still  remain;  the  wave  function  is  still  truncated  (now  at  z = -s 
instead  of  z = -s  + s^) , and  the  procedure  used  to  evaluate  the  shifts 
is  not  a standard  perturbation  technique. 

Rather  than  use  these  cutoff  methods,  we  may  employ  a consistent 
(though  still  approximate)  method  of  treating  the  image  problem  in 
metals.  The  simplest  such  technique  is  to  use  the  Thomas-Fermi 
model^ for  the  dielectric  response  of  the  metal.  In  this  theory,  the 
image  potential  is  determined  by  the  partial,  differential  equations 


V^(j)(r')  = 4iTq6(r')  - 4-rTq6(r'  - r)  , z'  > -s  , 
V^ct>(r')  = K^^4>(r’)  , z'  < -s  , 


(17) 


where 


2 _ 6TTNe^  _ r 12ff  1 

^ L(3.2)2/3“oJ 

with  N the  number  of  free  electrons  per  unit  volume  and  a the  Bohr 
radius.  The  total  electrostatic  energy  is  then  obtained  from 
equation  (17)  by  calculating  the  total  wor)c  necessary  to  bring  the 
charges  to  their  respective  positions  from  infinity. 


Pines,  Elementary  Excitations  in  Solids,  Benjamin,  New  York 
(1963),  96. 
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r 


5.2  Renormalization  Terms 


I 


t 

i 


1 


In  the  analysis  presented  by  Gadzuk,  first-order  perturbation 

theory  was  applied  only  to  the  interaction  potential  with  the  wave 

functions  normalized  in  the  region  z > s - s.  An  additional  shift 

c 

arises  if  one  considers  the  unperturbed  Hamiltonian 


H = T 
o 


(18) 


where 


T 


and  p is  the  momentum  operator,  and  if  one  evaluates  the  matrix  of  H 

o 

using  the  perturbed  wave  functions.  If  then  the  matrix  elements  of  H 

o 

for  large  s are  subtracted  out,  an  additional  shift  is  obtained,  or 


where 


and 


<s I H I s> 

AE  = 1 - <“|h^|“»  , 

n <s  s>  0 


oo 

/oo  on 

dz  If  i)i*H^ii.'dxdy 

~ Q ^ -00  — oo 


<S  S 


00 

/oo  no 

-// 


i|j*iJ/dxdy  . 


s -s 
c 


(19) 


(20) 


(21) 


By  using  the  wave  function  given  in  equation  (6) , the  result  is 

)]  / (22) 


AE  = 3.81a^  + 3.60a  ffafs  - s 

n L ’ c' 


1 
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where 


-2x  [j  x-x^+x^l 

-2x  li  3 x^  iJ' 
® (2'^4^'^2  "^2, 


Equation  (23)  was  used  to  calculate  Ae  for  the  two  cases  considered  by 

n , 

Gadzuk  (fig.  13) . A comparison  of  these  results  with  those  of  figure  4 

O 

shows  that  at  s = 3 A the  AE^  is  approximately  63  percent  of  the  result 
given  in  figure  4.  The  corresponding  quantity  at  s = 4 A is  still 
approximately  20  percent.  (Gadzu)c  chose  the  region  s = 3,  4 A,  for 
comparison  of  his  results  to  reported  data.)  In  this  region,  the 
contribution  of  Ae^^  is  certainly  not  negligible. 


Sj  DETERMINED  FRDM  S5r[(l+A)  l$  + 9^'  ' -J-3« 
14.4  2 


DISTANCE  OF  NUCLEUS  FROM  SURFACE  |A| 


Figure  13.  Energy  shift  due  to  wave-function  renormalization. 


6.  DISCUSSION  AND  CONCLUSION 

The  energy  shift  of  alkali  earth  atoms  near  metals  has  been 
calculated  based  on  a modified  version  of  Gadzuk's  theory.  This  energy 
shift  includes  image  effects  and  effects  due  to  the  ionic 
pseudopotential  and  the  state  of  the  ions  at  the  surface.  We  have  found 
several  errors  in  Gadzuk's  paper  which,  when  corrected,  alter  some  of 
the  consequences  of  that  work.  We  have  found  shifts  in  energy  of  about 
0.6  eV  at  distances  of  several  angstroms  from  the  metal  surfaces,  which 
is  a factor  of  two  higher  than  sliifts  found  by  Gadzuk.  Finally,  we  have 
shown  that  Gadzuk's  theory  is  incomplete  in  that  it  neglects 
wave-function  renormalization  effects,  and  we  have  proposed  an 
alternative  model,  the  Thomas-Fermi  model,  for  image  forces  in  the 
metal. 

What  is  of  primary  interest,  but  has  not  been  addressed  yet,  is  the 
transition  rate  of  the  electron  from  the  alkali  atom  to  the  metal.  To 
calculate  this  quantity  to  a sufficient  rate  or  precision,  the  potential 
seen  by  the  electron  must  be  known  quite  accurately.  We  conclude  that 
Gadzuk's  theory  does  not  permit  accurate  knowledge  of  the  potential  and 
that  a more  self-consistent  treatment  such  as  the  Thomas-Fermi  method 
must  be  used. 


LITERATURE  CITED 


(1)  J.  C.  Slater  and  K.  H.  Johnson,  Quantum  Chemistry  and  Catalysis, 
Physics  Today  (October  1974),  34. 

(2)  C.  A.  Morrison,  N.  Karayianis,  and  D.  E.  Wortman,  A Possible  Use  of 
the  Surface  States  of  Transition  and  Rare-Earth  Metal  Ions  in  the 
Theory  of  Catalysis,  Harry  Diamond  Laboratories  TR-1752  (April 
1976) . 


(3) 

C. 

A. 

Morrison 

, Phys . 

Lett . , 51A 

(1975) , 

49. 

(4) 

J. 

W. 

Gadzuk, 

Surface 

Science,  ^ 

(1967) , 

133. 

(5) 

J. 

w. 

Gadzuk , 

Surface 

Science,  6 

(1967) , 

159. 

(6)  F.  Herman  and  S.  Skillman,  Atomic  Structure  Calculations,  Prentice 
Hall,  Inc.,  Englewood  Cliffs,  HJ  (1963). 

(7)  M.  Abramowitz  and  I.  A.  Stegun,  ed. , Handbook  of  Mathematical 
Functions,  National  Bureau  of  Standards,  Gaithersburg,  MD  (1964). 

(8)  W.  A.  Harrison,  Pseudopotentials  in  the  Theory  of  Metals,  Benjamin, 
New  York  (1966)  . 

(9)  B.  G.  Wybourne,  Spectroscopic  Properties  of  Rare  Earths, 
Intersciences  Publishers,  New  York  (1965). 

(10)  F.  M.  Propst,  Phys.  Rev.,  129  (1962),  7. 

(11)  D.  Pines,  Elementary  Excitations  in  Solids,  Benjamin,  New  York 
(196  3)  , 96. 


SYMBOLS 


a 

a 

0 
B 

nm 

ft 

H 

m 

N 

P 

q 

Qi 

r 

R 

s 

s 

c 

V, 

1 

V 

o 
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lattice  constant 
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crystal  field  parameter 
Fermi  energy 
Planck's  constant 
single  electron  Hamiltonian 
electron  mass 

number  of  free  electrons  per  unit  volume 

momentum  operator 

electron  charge 

residual  charge 

ion  core 

ion  image 

distance  of  nucleus  from  surface 
cutoff  distance 
interaction  energy 

sum  of  Fermi  energy  and  work  function 
electrostatic  potential  energy 
surface-state  potential 
coordinates  normal  to  metal  surface 
Dirac  delta  function 
energy  shift 
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SYMBOLS  (CONT'D) 

energy  shift  due  to  renormalization 
surface-charge  density 
work  function 
complex  conjugation 
Laplacian  operator 
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APPENDIX  A. --ENERGY  SHIFT  IN  GADZUK'S  THEORY 


In  this  appendix,  we  derive  the  energy  shift  due  to  the  image 

charges  as  originally  presented  by  Gadzuk.'  We  set  up  a coordinate 

system  centered  about  the  ion  core  of  the  ad-atom,  with  the  z-axis 
normal  to  the  surface  of  the  metal.  In  this  coordinate  system,  the 

nucleus  is  at  (0,0,0)  (fig.  1 in  main  body  of  report),  and  the  electron 

is  at  point  r = (x,y,z)  . The  image  potential  ('^’imagg)  point  r due 
to  charge  q at  r'  with  an  ideal  metal  to  the  left  of  plane  z = -s  is 
given  by 


■ 

image 


(A-1) 


where 


e s, 
z 

(x’,y’,-z')  . 

If  we  have  an  assembly  of  charges,  the  total  potential  energy  of 
the  system  is  equal  to  the  total  wor)<  done  in  bringing  the  charges  to 
their  positions  from  infinity.  From  equation  (A-1),  we  calculate,  for 
the  system  of  figure  A-1, 


U. 

1 


(A-2) 


'j.  W.  Gadzuk,  Surface  Science,  ^ (1967),  159. 
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where 

is  the  total  potential  energy, 
q is  the  electronic  charge, 
r is  the  ion  core , 

s is  the  distance  of  the  nucleus  from  the  surface. 


In  equation  (A-2) , the  first  term  represents  the  direct  electron-ion 
interaction;  the  second  is  the  interaction  of  the  electron  with  its  own 
image;  the  third  is  the  interaction  of  the  electron  with  the  ion  image; 
and  the  last  is  the  ion- ion  image  interaction.  The  first  term, 
is  included  in  the  unperturbed  Hamiltonian,  while  the  last  term  does  not 
depend  on  the  electron  coordinates. 


Finally,  we  make  the  approximation  ]2s  + £|  - (2s  + z)  in  the  third 
term  in  equation  (A-2).  Therefore,  the  second  and  third  terms  in 
equation  (A-2)  are 


V. 

1 


4(s  + z) 


2s  + z 


(A-3) 


(where  is  the  interaction  energy)  which  is  the  form  used  by  Gadzuk. 

In  calculating  the  energy  shift  due  to  the  potential  of 

equation  (A~3) , Gadzuk  made  the  additional  assumption  that  the  electron 

wave  function  vanishes  for  z < s - s,  where  s , the  cutoff  parameter, 

c c 

is  given  by 

= -V  , (A-4) 

Z=S  -S  o 

c 

where  V is  the  sum  of  the  Fermi  energy,  E , and  the  work  function,  <}, . 
o r 

Equation  (A-4)  yields  the  following  result  for  s^ : 

Sc  = J |[(s  + X)  (s  + 9X)]  - s - 3xj  , (A-5) 


V.  (z) 
1 
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where 

X = qV<V  = 3.6A^  , 

o o 

V is  in  electron  volts, 
o 

A is  in  angstrom  units. 


The  electron  wave  function  is  taken  as 


i|)(r)  = N(1  - ar)  e 


-ar 


(A-6) 


vrtiere 

iji  is  the  wave  function, 

N is  the  normalization  constant, 

and  ij;  is  assumed  to  vanish  for  z < s - s.  Therefore,  N is  given  by 


oo  ou 

^.2,/  dz  / 


J /I  “*2ar 

pdp  (1  “ ar)  ‘•e 


(A-7) 


s -s 
c 


/ 2 

where  r = \^P  + z j . The  energy  shift  is  obtained  by  first-order  pert- 

urbation theory: 


oo  00 

AE^  = 2ttN^  J dzV^(z)  J pdpd  - ar)  , 


(A-8) 


®c“^ 


where  Ae^  is  the  energy  shift.  The  integral  over  p may  be  easily  per- 
formed by  making  the  substitution  p = r - z.  The  result  is 


AE.  = e2 


00 

dzF(|z|)  [25  V-"-  4(s  i"7)] 


(A-9) 


f 


s -s 
c 


dzF(  z ) 
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where 


. 1 -2az  /l 

F(z)  = e I ~ + az  - a 

4«2  '^2 


2z2  + 2 


(A- 10) 


and  the  z-integrals  may  be  performed  to  give 

e^a  1 

AE^  = |4g(2as)  (Ei(4as)  - Ei[2a(s  + s^)J) 


-4g (-2as) Ei (-4as)  - g (as) [e1 (2as)  - Ei(2as^)] 


+g  (-as)  Ei  (-2as)  + “ s^)  + 25(as)2 


-13a^ss  + 3(-_ 

c c 


as^)  2jexp  [-2a  (s  - 3^)^! 


(A-11) 


1 2 - expj^-2a(s  - ^ ^ ^ ^ 


s - s„) 


’]!■ 


where  g(x)  = ^4  ■*"  ^ “ ^2  + 2x^)e  This  is  the  form  used  in  the  main 

body  of  this  report.  In  equation  (A-11) , the  exponential  integral  is 


•OO 


Several  convenient  forms  for  computation  are  given  by  Abramowitz  and 
Stegun. 2 


2.M.  Abramowitz  and  I.  A.  Stegun,  ed. 
Functions , National  Bureau  of  Standards , 


, Handbook 
Gai thersburg , 


of  Mathematical 
MD  (1964)  . 
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s 

c 

U. 

1 

V, 

1 

V 

o 

Ae. 


image 

4> 
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unit  vector  in  z direction 
Fermi  energy 
electronic  charge 
ion  core 

distance  to  point 

distance  of  nucleus  from  surface 

cutoff  parameters 

total  potential  energy 

interaction  energy 

sum  of  Fermi  energy  and  work  function 

energy  shift 

image  potential 

work  function 

wave  function 
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APPENDIX  B.— POTENTIAL  DUE  TO  IONS  IN  SOLID 


B-1,  INTRODUCTION 

In  this  appendix,  we  discuss  two  contributions  to  the  electrostatic 
potential  outside  a metal  surface.  One  contribution  is  due  to  the 
lattice  of  point  charges  and  the  uniform  background  of  the  electron  gas. 
The  second  contribution  is  due  to  the  asphericity  of  the  ions 
constituting  the  solid;  this  asphericity  depends  on  the  electronic  state 
of  the  ions . 

B-2.  POINT-CHARGE  LATTICE  -SUM 

In  discussing  the  first  contribution  to  the  potential,  we  consider 
a rectangular  lattice  with  lattice  constants  a and  b.  It  is  assumed 
that  a point  charge,  Qq,  is  situated  at  the  center  of  each  unit  cell. 
At  the  same  time,  we  consider  a uniform  background  of 
magnitude  o = -Qq/ab.  It  is  necessary  to  consider  both  these 
contributions  together  in  order  to  obtain  finite  results. 

We  consider  only  the  first  layer  of  point  charges  in  the  metal.  We 
may  do  so  since  screening  by  the  electron  gas  is  expected  to  reduce  the 
contributions  of  other  layers  considereibly , and  even  without  this 
screening  the  contributions  due  to  the  other  layers  are  small,  as  shown 
by  detailed  calculations.  We  assume  that  the  first  layer  is  not 
screened  at  all. 

To  properly  handle  the  divergence  problems  inherent  in  the 
calculation  of  the  potential,  we  consider  a finite  system  consisting  of 
(2Ri  + 1)  unit  cells  in  the  x-direction  and  (2R2  + 1)  unit  cells  in  the 
y-direction.  The  origin  of  coordinates  is  taken  at  the  center  of  the 
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lattice;  we  assume  one  ion  at  the  center  of  each  unit  cell.  The  total 
electrostatic  potential  energy  (V^)  of  a charge  -q  at  the  point  (x,y,z) 


is  then  given  by 


+Rj  '*■^1 


V (x,y,z)  = ^ Y.  7 ^ ^ ; 77 

P P=-Rl  q'=-R2  - pa)^  + (y  - q'b)^  + z^J 

a (Rj  +^5)  b (R2+*5) 

+ dx'  dy' 

-a(Ri+S)'  -b(R2+*5) 


tT  • VB-j.; 

(X  - x')^  + (y  - y')^  + 

The  first  term  in  equation  (B-1)  is  due  to  the  point-charge  lattice,  and 
the  second  is  due  to  the  uniform  background.  We  are  interested  in 
evaluating  equation  (B-1)  in  the  limit  Ri,R2  “*■ 

Equation  (B-1)  may  be  transformed  to  a more  convenient  form  by 
making  the  substitution 


pa)^  + (y  - q'b)^  + ^ 


oo  . 

^ J t dt  expj-t[(x  - pa)  ^ + (y  - q'b)^  + z^]j  . 


Now  we  must  consider  sums  of  the  form 


S(t)  = Y exp[-t(x  - pa)^] 
p=-Rj 
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and  use  the  method  of  converting  such  sums  given  by  Whittaker  and 
Watson.'  Consider  the  function  f(z)  = exp  [-t(x  - za)^J  in  the  complex 
plane.  By  considering  the  contour  shown  in  figure  B-1,  we  may  show  that 


:(t)  = f 


r dzf (z) 

f 2iriz  , 
C e - 1 


lim  r dx'f  (X')  _ f dx 

e-K)  J 2Tii(x'-ie)  ^ J 2Tri(x 

■1  . e - 1 "A  1 6 


dx'f (x' ) 
(x*  +ie) 


-R^-h 


-Ri-S 


Now  we  may  expand  the  denominators  in  equation  (B-4)  and  take  the  limit. 


We  obtain 


R1+*5 

OU 

(t)  = L J.  dx'f(x')  e 

p=-oj  -Rj-% 


-2iripx' 


Now  we  consider  the  sum  in  equation  (B-1) . We  substitute 
equation  (B-2)  into  equation  (B-1)  and  use  equation  (B-5)  to  obtain 


V^(x,y,z)  = - 


00  — V 


dt  , o , 
: — exp (“tz^ ) 


p,q*=-oo  o 


J f dy'e-2^"^'y' 


-Ri->5 


-Rt-^j 


• expj-t  [(x  - x'a)^  + (y  - y'b)^]| 


a(Rj+Jj)  b(R2+*5) 

I / 

— a(R|+*3)  ~b(R2+^) 


[(x  - x')^  + (y  - y')^  + z^]*’ 


'e.  T.  Whittaker  and  G.  N.  Watson,  A Course  of  Modern  Analyses, 
University  Press,  Cambridge,  U.K.  (1950),  124. 
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Figure  B-1.  Contour  used  for  integral. 


Now  we  consider  this  expression  in  the  limit  as  Rj  and  R2  become 
infinite.  The  integrals  over  t converge  in  the  limit  except  for  the 
integral  with  p = q'  = 0.  We  consider  this  term  separately.  We  have 


^00  — 

Qq2  r 

V7 

exp (-tz2) 

R1+>3 

/ ax' 

-Rl-', 

R2+‘5 

/ ay 

-R2-*5 

■ exp|-t  [(x 

- x' a) 2 + 

(y  - y’b)2]j= 

-0q2  J 

dx-  J 

dy ' 

1 

-Rl-h 

-R2-"5 

r,x 

- x'  a)  2 + (y 

and  this  contribution  cancels  the  last  integral  in  equation  (B-6) . In 
the  remainder  of  equation  (B-6) , we  may  immediately  take  the  limit  Ri , 
R2  00.  have 


, -2Ttipx'  — 
dx'e  ^ exp  -t 


(X  - x-a)2]  = i^expf-  ^ - 2liEx]  , 

a II  t 2 S 

L ta^  J 


(B-8) 


with  a similar  result  for  the  y'  integration.  Therefore,  we  obtain 


Vp(x,y,z)  . - g- 


exp 


^f\- 


3/2 


dt  exp 


(B-9) 
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where  the  prime  sign  on  the  summation  in  equation  (B-9)  indicates  that 
the  term  p=q'=0  is  to  be  omitted.  The  integral  in  equation  (B-9)  may 
be  performed  to  obtain 


V^(x,y,z) 


ab  L—J 

p,q'=_oo 


exp 


-2.i(Ei  + ^)  _ 27t|z| 


(B-10) 


This  is  the  final  result. 


B-3.  CONTRIBUTION  OF  SURFACE  STATES 


We  now  consider  those  contributions  to  the  surface  potential  which 
depend  on  the  states  of  the  ions  at  the  surface  of  the  metal.  Consider 
an  ion  whose  center  is  on  the  lattice  at  the  point 

~i  " ®x^^  ' (B-11) 

where  R is  the  vector  from  the  origin  to  the  ion,  e and  e are  unit 

X y 

vectors  in  the  x and  y directions,  respectively,  and  an  electron 

associated  with  the  ion  at  distance  r'  relative  to  R . Then  the 

~ ~i 

potential  at  point  r due  to  this  electron  is  given  by 


1 


q 

IR.  + r'  - r 

:~i  ~ ~ 


-E 


q(r')  (-1) 


(r’)C, 


(B-12) 


where  f and  R.  - r are  unit  vectors  in  the  direction  of  r'  and  R - r 
respectively,  and  the  C^,^  are  related  to  the  spherical  harmonics  by 


41 


APPENDIX  B 


c =r4TT/(2Jl  + 1)1^^  Y as  defined  elsewhere.^  Now  we  assume  that  the 
Hm  ■*  dm 

ion  is  in  a certain  state  and  that  the  potential  due  to  the  ion  may  be 
obtained  by  summing  equation  (B-12)  over  all  the  electrons  in  the  ion 
and  taking  the  expectation  value  of  the  result  over  the  ion  state.  We 
therefore  obtain,  for  the  potential  at  r, 


i.  (r) 


-E 

d,m 


(-1)  q+  C„  (R. 
dm  dm'~i 


R.  - r 


d+1 


(B-13) 


where 


q„  = - ^ 

dm  e 


<’Jj|r'  Cj^^(r’ ) 


Quantity  q„  is  thus  a function  of  the  state  of  the  surface  ion. 
The  d = 0 part  of  equation  (B-13)  is  the  point-charge  lattice  sum 
considered  in  section  B-2. 

For  the  other  terms  in  equation  (B-13) , we  assume  that  all  the 

surface  ions  are  in  the  seune  state.  We  therefore  obtain,  for  the 

potential  energy,  V , of  an  electron  at  r, 
s ~ 


The  sum  over  i in  equation  (B-14)  is  a lattice  sum  similar  to  that 
considered  earlier.  It  may  be  converted  to  a more  convenient  form  by 
the  same  techniques  illustrated  in  section  B-2  for  the  point-charge  sum. 
The  result  is 


^C.  A.  Morrison,  N.  Karayianis , and  D.  E.  Wortman,  A Possible  Use  of 
the  Surface  States  of  Transition  and  Rare-Earth  Metal  Ions  in  the  Theory 
of  Catalysis,  Harry  Diamond  Laboratories  TR-17S2  (April  1976). 
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SYMBOLS 


a,b 


lattice  constants 


e , e unit  vectors  in  x and  y directions 

X y 


R. 


V 

P 


V 

s 

x,y,z 
x' ,y* ,z' 


momentum  operator 
point  charge 
point 

unit  ve'tnr  in  direction  of  r' 
unit  vector  in  direction  of  - r 
vector  from  origin  to  ion 
electrostatic  potential  energy 
potential  energy 

coordinates  normal  to  metal  surface 


o 

<(>. 

1 

<}> 


magnitude 
potential 
work  function 
wave  function 
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APPENDIX  C.  ENERGY  SHIFT  DUE  TO  POINT-CHARGE  AND  SURFACE-STATE 
POTENTIALS 

In  calculating  the  energy  shift  due  to  the  ions  in  the  metal,  for 
each  type  of  ion  on  the  surface,  we  must  consider  sums  of  the  type  shown 
in  equations  (B-10)  and  (B-15)  in  appendix  B of  this  report.  The 
coordinates  of  the  ion  in  the  unit  cell  relative  to  some  arbitrary 
origin  are  contained  in  the  x and  y in  these  equations.  There  are  three 
contributions  to  the  x and  y: 

a.  Coordinates  of  the  nucleus  of  the  ad-atom  relative  to  the 

origin  of  the  lattice,  x , y 

o o 

b.  Coordinates  of  the  lattice  ions  in  the  unit  cell  relative  to 
the  origin,  x^,  y^ 

c.  Coordinates  of  the  ad-atom  electron  relative  to  the  nucleus, 

X,  y 


Thus,  the  term  [(px/a)  + (q'y/b)]  in  equations  (B-10)  and  (B-15)  is 
to  lae  replaced  by 


E 

a 


(x  + x - x.)+^(y  + y 


For  the  PrSb  lattice  in  the  [100]  plane,  the  ad-atom  is  directly 
above  a Pr  ion.  In  this  example,  the  lattice  is  square  (a  = b) ; for 
the  [100]  plane,  Pr  ions  are  at 


xi  = 0,  yi  = 0 , 

X2  = a/2  , y2  = a/2  . 


(C-1) 


PAlSiBLAliK-iyCT  FIL’/SD 


47 


APPENDIX  C 


The  Sb  ions  are  at 


X3  = 0,  y3  = a/2  , 
X4  = a/2,  yi*  = 0 . 


(C-2) 


Since  the  ad-atom  is  directly  over  a Pr  ion,  x = y =0.  There  is  a 

o o 

sum  over  i (1,  2,  3,  4)  in  equations  (B-10)  and  (B-15)  ; the  q are 

urn 

functions  of  i.  For  PrSb,  we  take  q„  (1)  = q (2)  and 

Hm  {.m 

~ ignore  the  surface-state 

contribution  and  consider  simply  the  point-charge  potential.  We  have  Pt 
ions  at  the  [100]  plane: 


XI  = 0,  yj  = 0 , 


X2  = a/2,  y2  = a/2 


(C-3) 


and  for  W we  have  one  ion  per  unit  cell  for  the  [100]  plane,  at 

xi  = yi  = 0. 

For  these  simple  lattices,  the  sum  over  i can  be  simplified.  For 
Pt,  we  take  only  terms  with  (p  + q')  even  and  (p  + q')  odd  in 
equations  (B-10)  and  (B-15)  and  multiply  the  result  by  2.  For  PrSb,  we 
take  the  sum  of  the  (p  + q')  even  and  (p  + q')  odd  terms  for  the  Pr  ions 
and  the  difference  of  the  (p  + q')  even  and  (p  + q')  odd  terms  for  the  Sb 
ions  and  multiply  both  Pr  and  Sb  by  2.  Doing  so  removes  the  i sum  and 
eliminates  the  dependence  of  equations  (B-10)  and  (B-15)  on  the 
unit-cell  coordinates. 

To  obtain  the  -'energy  shift  due  to  these 
ejqjectation  value  of  thQ  function 

g(r)  = exp  |-2Tii  - 


terms,  we  consider  the 
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when  z is  now  nveasured  relative  to  the  nucleus  of  the  ad-atom.  We  need 
the  integral 

J j dz  J pdp  (1  - ar)2e“^°‘'^g(r)  i — j-G(s  - s^,  a,p,q').  (c-5) 

0 S “S  0 

c 


The  angular  integral  is 


2n 


J~  d<P  exp  1^- 


2n 


EX 

a 


2nJ 


(C-6) 


where  is  the  Bessel  function  of  zero  order  and  the  p integrals  may  be 
performed  by  substituting  p = (r^  - z^)^  into  equation  (C-5).  We  use 
the  relations 


APPENDIX  C 


r 


Thus,  substituting  equations  (C-7)  to  (C-9)  into  equation  (C-5) , we  are 
left  with  the  z-integral,  which  is  an  integral  of  an  exponential  in  z 
multiplied  by  a polynomial  in  z.  The  result  is 

G(s  - s^,a,p,q')  = 4a3  ^ [4(R  - Y)^A  + 2a(R  - - ) + 

l8(R  - 


5 ~ S ^ 

2«2(r  - ,)c  + 3a^o]  + ^ [2a(R  - 0^8  + 2a2(R  - ,)C  + Ba^oJ  (c-10) 

4 (R  - y)^ 
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We  may  use  the  above  to  calculate  the  energy  shifts.  Defining 


- G(s  - s^,a,p,q')Jy^ 


2 - e 


■2a  (s-s  ) I , 3 , , , , , 1 

c 1^1  -t-  -a(s  - Sc)  + a2(s  - s^)  ^ (s  - Sc)Ht  , 


we  have  for  the  point-charge  potential  shift,  AE  , 

P 


AE  =-2^ 
p ab 


H(s  - s^,a,p,q' ) 


exp  <-27Ti 


cind  for  the  surface-state  potential  shift,  AE  , 

s 


AE  = - 3-  Y.  /^+ 

s ab  . ^ ^Jlm  2 v-2  / 


m -I 


H(s  - s^,  a,  p,  q') 


I Tp(Xc  + Xi)  q'(yo+yi)l  /2  ,2\^| 

exp  \-27Ti  + - 2tt  z /—  + 1 , 

( L » K J I \ 


, m < 0 . 


(C-14) 
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SYMBOLS 


lattice  constants 


Bessel  function  of  zero  order 


momentum  operator 
electronic  charge 
charge  on  lattice  ion 
distance  of  nucleus  from  surface 
cutoff  distance 

coordinates  of  ad-atom  electron  relative  to  nucleus 
coordinates  of  lattice  ions  in  unit  cell  relative  to  origin 
coordinates  of  nucleus  of  ad-atom  relative  to  origin 
point-charge  potential  shift 
surface-state  potential  shift 
potential 
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